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ABSTRACT 

We have used the Karl G. Jansky Very Large Array (VLA) to carry out a multi-epoch radio continuum 
monitoring of the Orion Nebula Cluster (ONC) and the background Orion Molecular Cloud (OMC) 

(3 epochs at Q-band and 11 epochs at Ka-band). Our new observations reveal the presence of 19 
radio sources, mainly concentrated in the Trapezium Cluster and the Orion Hot Core (OHC) regions. 

With the exception of the Becklin-Neugebauer (BN) object and the source C (which we identify here 
as dust emission associated with a proplyd) the sources all show radio variability between the different 
epochs. We have found tentative evidence of variability in the emission from the massive object related 
with source 1. Our observations also confirm radio flux density variations of a factor >2 on timescales 
of hours to days in 5 sources. One of these flaring sources, OHC-E, has been detected for the first 
time. We conclude that the radio emission can be attributed to two different components: i) highly- 
variable (flaring) non-thermal radio gyrosynchrotron emission produced by electrons accelerated in the 
magnetospheres of pre-main sequence low-mass stars; ii) thermal emission due to free-free radiation 
from ionized gas and/or heated dust around embedded massive objects and proplyds. Combining our 
sample with other radio monitoring at 8.3 GHz and the X-ray catalog provided by Chandra, we have 
studied the properties of the entire sample of radio/X-ray stars in the ONC/OMC region (51 sources). 

We have found several hints of a relation between the X-ray activity and the mechanisms responsible 
for (at least some fraction of) the radio emission. We have estimated a radio flaring rate of ^0.14 flares 
day“^ in the dense stellar cluster embedded in the OHC region. This suggests that radio flares are 
more common events during the first stages of stellar evolution than previously thought. The advent 
of improved sensitivity with the new VLA and ALMA will dramatically increase the number of stars 
in young clusters detected at radio wavelengths, which will help us to improve our understanding of 
the origin and nature of the radio emission. 

Subject headings: 


1. INTRODUCTION 

High energy processes during the first evolutionary 
stages of star form ation are responsible for both radio 
and X-ray emission (jFeigelson fc Montmerlelfl999[ ) . Low- 
mass pre-main sequence (PMS) stars are well known 
strong X-ray emitters. Their enhanced magnetic activ¬ 
ity with respect to more evolved stars produces violent 
reconnection events in the corona of the stars, where the 
plasma heated to high temperatures strongly emits vari¬ 
able X-ray emission. Massive stars also emit X-ray radi¬ 
ation, usually related to wind shocks. 

Our understanding of the X-ray emission from young 
stars has dramatically increased in the recent years 
due to Chandra an d XMM-Newton (|Getman et al.l[2QQM 
lArzner et all I2QQ7I ) . X-ray observations have revealed 
thousands of PMS stars in tens of stellar clusters, result- 
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ing in good constraints on their X-ray properties such as 
plasma temperatures, levels C )f variability, lum inosities 
and X-ray flare rate (see, e.g.. lWolk et al.l[2QQ5[ ). 

In contrast, the physics associated with the radio 
events (nature and origin of the emission, variability, 
timescales, flaring rat e) are still poorly constrained. 
iDrake fc Linskvl (|l989f ) proposed that radio flares might 
be produced by the same coronal activity that is re¬ 
spon sible for bright X-ray emission (see review bv lGiidel 
l2QQ2f ) . It would be expected then that electrons spiraling 
in the magnetic field of the corona produce non-thermal 
and highly variable gyrosynchrotron radiation. More¬ 
over, ionized material in the vicinity of stars, in circum- 
stellar disks or envelopes or at the base of bipolar out¬ 
flows, also produce thermal free-free (bremsstrahlung) ra¬ 
diation. 

Long-term radio variability on timescales of months 
to years has be e n observed in star- f orming regions 
(iFelji et al.l Il99^ IZapata et al.l I2QQ4 iForbrich et al.l 
l2QQ7l : IChoi et all l2QQ8f ). However, it is still not clear 
whether these variations are caused by long-term mech¬ 
anisms, or they are indeed produced by a sequence of 
events occurring on shorter timescales. Systematic ob¬ 
servations looking for short-term variability are required 
to answer t his questio n. _ 

Recently, iLiu et aP (j2QI4l ) detected radio variability on 
hour timescales in the young stellar cluster R Coronae 
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Table 1 

Multi-epoch VLA radio continuum observations. 


Config. 

Project 

Band 

Freq. 

BW 

Epoch 

JD 

Pointing center 

Obs. 

beam 

RMS 

Primary 

Gain calibrator 


ID 

name 

(GHz) (MHz) 



RAj2ooo DECj2ooo length 

(" X ") 

(mJy) 

fiux density fiux density (Jy) 








5h 35m 

-5° 22' 

(min) 



calibrator 

(J0541-0541) 

B 

AJ356 

Q 

45.6 

25 

2009 Mar 9 

2454900 

14.60s 

30.0" 

120 

0.23x0.15 

0.45 

J0137+3309 

0.63 

B 

AJ356 

Q 

45.6 

25 

2009 Mar 19 

2454910 

14.60s 

30.0" 

240 

0.22x0.15 

0.37 

J0137+3309 

0.63 

D 

AR712 

Q 

43.3 

100 

2009 Dec 22 

2455188 

14.50s 

31.0" 

60 

1.9x1.4 

0.54 

J0137+3309 

0.46 

C 

lOB-175 

Ka 

33.6 

256 

2010 Oct 24 

2455494 

14.50s 

30.0" 

30 

1.3x0.65 

0.36 

J0542+4951 

0.68 

G 

lOB-175 

Ka 

33.6 

256 

2010 Nov 23 

2455523.76 

14.50s 

30.0" 

30 

0.92x0.60 

0.40 

J0542+4951 

0.68 

C 

lOB-175 

Ka 

33.6 

256 

2010 Nov 23 

2455523.89 

14.50s 

30.0" 

30 

0.75x0.56 

0.44 

J0137+3309 

0.69 

G 

lOB-175 

Ka 

33.6 

256 

2011 Jan 8 

2455569 

14.50s 

30.0" 

30 

0.86x0.56 

0.39 

J0137+3309 

0.62 

CnB-B 

lOB-175 

Ka 

33.6 

256 

2011 Feb 8 

2455601 

14.50s 

30.0" 

30 

0.30x0.24 

0.32 

J0137+3309 

0.58 

B 

lOB-175 

Ka 

33.6 

256 

2011 Mar 28 

2455649 

14.50s 

30.0" 

30 

0.28x0.22 

0.22 

J0542+4951 

0.50 

BnA 

lOB-175 

Ka 

33.6 

256 

2011 May 27 

2455709 

14.50s 

30.0" 

30 

0.62x0.077 

0.22 

J0542+4951 

0.54 

BnA-A 

lOB-175 

Ka 

37.5 

128 

2011 Jun 4 

2455717 

14.50s 

30.0" 

120 

0.096x0.073 

0.13 

J0542+4951 

0.50 

BnA-A 

lOB-175 

Ka 

30.5 

128 

2011 Jun 4 

2455717 

14.50s 

30.0" 

120 

0.13x0.083 

0.10 

J0542+4951 

0.52 

A 

lOB-175 

Ka 

33.6 

256 

2011 Jun 11 

2455724 

14.50s 

30.0" 

30 

0.14x0.062 

0.16 

J0137+3309 

0.54 

A 

lOB-175 

Ka 

33.6 

256 

2011 Jul 09 

2455752 

14.50s 

30.0" 

30 

0.091x0.060 

0.14 

J0542+4951 

0.48 
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Fig. 1.— Position of the 19 radio sources detected in our monitoring program, overplotted on the R-band ACS/WFC HST image (left 
panel) and on the K-band 2 Micron All Sky Survey (2MASS) image (right panel). The dotted circles indicate the regions of the field where 
the primary beam responses are >0.15 at Q-band (43.3 GHz) and Ka-band (33.6 GHz). 


Australis. The most powerful radio flares^ have been 
serendipitously reported so far toward the Orion Neb¬ 
ula Cluster (P NC) and the backgr ound Orion Molecular 
Cloud (OMCb [Bower et all ()2QQ3[ ) reported a strong ra¬ 
dio flare at 86 GHz arisin g from a PMS star in the ONC, 
and iForbrich et aTI (|2QQ8[) presented an even stronger ra¬ 
dio flare at 22 GHz, originating from a young star deeply 
embedded in the OMC previously detected through its 
X-ray emission. The low number of observed event s could 
indic ate that radio flares are a rare phenomenon (|Andrel 
Il996[ ) , but at the same time prevents a proper statistical 
analysis of short-term variability phenomena. 

Since the typical timescales of radio variability are 

® In this work we will use the term flare to refer to fiux density 
variations of a factor of >2 on timescales from hours to days. 


poorly known, we have carried out a monitoring pro¬ 
gram comprising various cadences ranging from 3 hours 
to several months. The new capabilities of the Karl G. 
Jansky VLA now allows the scheduling of multiple, short 
snapshots with good sensitivity in a reasonable amount 
of observing time. The only two examples of powerful 
radio flares are located in Orion, and this region also har - 
bors a rich cluster of low-mass stars (|Rivilla et al.ll2Ql§) : 
hence, this region is an excellent target for the detection 
of many sources in a single pointing. 

We have carried out a multi-epoch radio continuum 
monitoring of the ONG/OMG region using the Karl G. 
Jansky Very Large Array (VLA). This is the first radio 
monitoring at high centimeter frequencies in Orion, with 
3 epochs at Q-band and 11 epochs at Ka-band. Our data 
allow us to study for the first time both the short (hours 
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to days) and long (months) timescale variability of the 
radio sources in Orion. 

The paper is laid out as follows. In Section [2] we present 
the details of the observations. In Section [3] we show the 
results of the monitoring at Ka and Q bands. In Section 
[U we compare our results with a previous monitoring at 
lower frequency (8.3 GHz). We also compile the full sam¬ 
ple of radio/X-ray sources in the ONC/OMC region, and 
then compare the radio and X-ray properties with the 
aim of better understanding the link between radio and 
X-ray emission. In Section [5] we discuss in more detail 
the new radio flaring source detected by our monitoring, 
OHC-E. In Section [6] we analyze in particular the radio 
variability observed in the binary system 0^ Ori A. In 
Section [3 we estimate the rate of radio flaring activity 
of young stars in the Orion Hot Core (OHC). Finally, in 
Section [8] we summarize the conclusions of our work. 

2. OBSERVATIONS AND DATA REDUCTION 

The radio observations of the ONC/OMC region were 
made with the VLA in the A, B, C and D configura¬ 
tions at Q-band (3 epochs) and Ka-band (11 epochs) 
from 2009 to 2011. Table [T] summarizes the observa¬ 
tional details for each of the epochs: array configuration, 
project ID, band name, frequency, observing bandwidth, 
date, pointing center, observation length (including all 
calibrations), synthesized beam, RMS noise at the cen¬ 
ter of the primary beam, primary flux density calibrator 
used, and the flux density of the complex gain calibra¬ 
tor. The separation between the epochs is different: the 
shortest one is only 3 hours, others are separated by sev¬ 
eral days, and some other by several months. This will 
allow to trace radio variability at different timescales. 
The first two observations at Q-band (45.6 GHz) were 
made in spectral line mode using the old VLA correlator, 
single polarization. The third Q-band observation (45.6 
GHz, AR712) used the standard dual polarization contin¬ 
uum set-up with the old correlator (100 MHz total band¬ 
width). For the data taken under observing code 1 OB-175 
the new WIDAR correlator was used, with two 128 MHz 
sub-bands centered on 33.56 GHz placed contiguously in 
frequency for all observations except those of 2011 Jun 
4, for which the sub-bands were separated by 7 GHz as 
noted in Table [H The observations at Q-band were re¬ 
duced using the Astronomical Image Processing System 
(AIPS) package. The observations at Ka-band were re¬ 
duced using the VLA Calibration Pipeline^, which uses 
the Common Astronomy Software Applications (CASA) 
package^^. 

Phase self-calibration was performed where possible. 
However, the tropospheric phase stability of the obser¬ 
vations on 2011 Jun 4 was poor, and the data were 
of insufficient signal-to-noise ratio to enable phase self- 
calibration on short enough timescales to correct for the 
resulting decorrelation. While the typical uncertainty in 
the absolute flux density scale is 10% at these frequencies, 
the uncertainty in the absolute flux density scale for the 
2011 Jun 4 data is increased to 20%. Images were made 
using CASA, applying an inner uv-cutoff (>50 kA) to fil¬ 
ter the extended emission from the foreground HII region 

® https://science.nrao.edu/facilities/vla/data- 

processing/pipeline 

http://casa.nrao.edu 


ionized by the Trapezium cluster, and natural weighting 
to maximize sensitivity. The images are corrected for 
the response of the primary beam before being used for 
photometry. The FITS files of the reduced images are 
available in the electronic version of the paper. 

3. RESULTS: RADIO STELLAR POPULATION DETECTED 
AT HIGH CENTIMETER FREQUENCIES 

The large fields of view (FoVs) covered by the images 
made here, especially in the VLA’s B and A configura¬ 
tions, requires a rigorous criterion for the detection of 
sources. Assuming a gaussian noise, and approximat¬ 
ing the number of potential sources as the ratio between 
the area of the FoV and the solid angle of the beam 

(7r^minor^major/4I/n2), WC WOUld CXpCCt < 0.3 SOUrCCS 

above 5 times the local RMS noise (5cr) in the worst case. 
Therefore, we only report detections with flux densities 
> 5(7. 

We detected a total of 19 sources (Table[2]), 18 of which 
have been previously detected by the radio monitorings 
at l ower frequenc i es (|Felli et al.l ll99^ at 5 and 15 GHz 
and IZapata et al.ll2QQ4l at 8.3 GHz). Our observations 
have revealed the presence of a new radio source, here¬ 
after OHC-E, detected in the OHC region in two different 
epochs. In Fig.[T]we show the positions of all the detected 
sources, overplotted on the R-band Advanced Camera for 
Surveys/Wide Field Channel (ACS/WFC) Hubble Space 
Telescope (HST) image and on the infrared K-band im¬ 
age from the 2 Micron All Sky Survey (2MASS). The 
sources are mainly concentrated in the OHC region and 
the Trapezium Cluster, which h arbors the two hig hest 
stellar densities within our FoV (|Rivilla et al.ll2Ql^ . 

To measure the flux densities of the sources, we use 
the AIPS task JMFIT^^. We also add an absolute uncer¬ 
tainty of 10% in quadrature (20% for the 2011 Jun 4 ob¬ 
servations due to the poorer phase stability of those data; 
see Section HD- Table El summarizes the results. Only 
the sources BN and source I are detected in all epochs, 
showing nearly constant flux densities, and verifying the 
reproduceability of the flux density scale. Source C is 
consistent with a constant flux density, although is not 
detected in all epochs. The other sources exhibit clear 
flux density variation between epochs. In the case of 
non-detections, we quote 3(7 upper limits. 

3.1. Long-term variability: month timescales 

In this section we study the behavior of the radio emis¬ 
sion throughout the full monitoring. In Figs. [3] and |4] we 
show the measured integrated flux densities of all sources 
during the different epochs of our monitoring for the Q- 
band and Ka-band, respectively. We note that in the fig¬ 
ures the observation on JD 2455188 (2009 Dec 22) at Q- 
band corresponds to a slightly different frequency (43.3 
GHz) from the first two epochs (45.6 GHz). Also, the 
flux density from the observation on JD 2455717 (2011 
June 4) at Ka-band (shown in Fig.|4]) corresponds to 30.5 
GHz, and not to 33.6 GHz as for the others. 

With the aim of quantify the radio variability we study 
two parameters: i) the standard deviation AF of the 
flux densities from the average flux Fav, which measure 
the absolute variability; and ii) yd, which is defined as 

JMFIT fits a 2 dimensional gaussian models to the sources by 
least-squares. 
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Fig. 2.— Images of the 19 radio sources detected in our radio monitoring. The name of the source is indicated in the top left corner of 
each panel. In the lower right corners we indicate the epoch of the observation (we have selected the epochs with the best spatial resolution 
and good detections, and in some cases we show the detections in two different epochs with different colors). The first contour level is 3(7 
and the step between successive contours is 2(7, with the exception of BN, I, 25 and 12 (which have steps of lOcr), and the 2009 March 19 
epoch for OHC-E (steps of her). The open blue squares indicate the position of known radio sources detected bv IZaoata et al.l (120041) at 
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Fig. 3.— Q-band light curves for those sources that fall within 
the Q-band primary beam. In the case of non-detections, 3a upper 
limits are indicated with triangles. The flux densities have been 
normalized by the average value (Tav) calculated with the positive 
detections (or the average value of the upper limits if the source is 
not detected). The error bars indicate the flux density uncertain¬ 
ties. The red horizontal line indicates the value at which the flux 
density is equal to the average. The first two observations were 
carried out at 45.6 GHz (black filled symbols), while the last one 
was carried out at 43.3 GHz (black open symbols). For source G, 


Table 2 

Positions of the radio sources detected in our 
Ka-band and Q-band monitoring. 


Source 

P^J2000 

5 35 

DEC J2QQQ 
-5 

BN 

I4.II 

22 22.69 

I 

14.51 

22 30.58 

OHG-E 

14.73 

22 29.83 

D 

14.90 

22 25.38 

na 

14.35 

22 32.89 

H 

14.50 

22 38.76 

A 

11.80 

21 49.29 

G 

14.16 

23 01.04 

F 

18.37 

22 37.43 

G 

17.95 

22 45.42 

E 

16.96 

22 48.78 

15 

16.07 

23 07.12 

6 

16.75 

23 16.44 

7 

16.28 

23 16.58 

25 

15.77 

23 09.86 

12 

15.82 

23 14.00 

II 

15.84 

23 22.40 

16 

16.33 

23 22.54 

5 

16.85 

23 26.31 


^ This source was ca lled L bv IGar ayj_([jj 87lb but its more 
common name is n (IMenten &: ReidllI995lb 


P = AF/Fav, following iFelli et abl (jl993f ). which mea¬ 
sures the relative variability. Given that some sources 
remain undetected in some epochs, and hence have flux 
densities below the sensitivity limit, we consider in the 
calculation of AF (and hence of P) the positive detec¬ 
tions as well as the lowest upper limit. In Table H] we 
show the values of Fav, AF and P for the sources de¬ 
tected. 

The sources BN and I are associated with massive 
stars (|Ruid et al.l I2QQ7L iGoddi et al.ll2Qll[ ). and are ex¬ 
pected to emit mainly thermal (and constant) emission 
arising from ionized gas surrounding the central object. 
Our monitoring at Q-band and Ka-band (Figs. [3]and|4]) 
shows indeed that the flux density of BN is nearly con¬ 
stant, with a very low month-timescale radio variability 
parameter at 33.6 GHz of /3=0.06. 

In the case of source I the radio variability is higher, 
/3=0.13. IZapata et al.l (|2QQ4f ) also reported some vari¬ 
ability towards this source, with flux density varia tions 
of a factor ^2. Furthermore. iPlambeck et ~ahl (|2Q13[ ). us¬ 
ing observations separated by 15 yr also found evidence 
of a gradual flux density increase from source I with re¬ 
spect the more steady flux density of BN (see their Fig. 
5). From our Ka-band monitoring, we have studied the 
month-timescale evolution of the ratio between the flux 
densities of sources I and BN. Fig.[5]shows that this ratio 
exhibits variations larger than the statistical uncertain¬ 
ties (we do not include the uncertainty in the absolute 
flux density scale in calculating the error in these ratios), 
suggesting that real variability is present. The origin of 
this confirmed long-term variability toward source I could 
be due to ion ization of infalling accretion flows onto the 
mass ive star (jGalvan-Madrid et al.l I2Q11I: iDe Free et al.l 

l2fnl) . 

Source G is only detected in the epochs for which the 
VLA was in its most compact configurations. This sug- 
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Fig. 4. — Ka-band light curves of all the 19 sources detected throughout the Ka-band monitoring. In the case of non-detections, 3a 
upper limits are indicated with triangles. The flux densities have been normalized by the average value (Tav) calculated with the positive 
detections (or the average value of the upper limits if the source is not detected). The error bars indicate the flux density uncertainties. 
The red horizontal line indicates the value at which the flux density is equal to the average flux density. The dotted line merely joins 
the flux densities at different epochs, and is not indicative of the evolution of the ffux density between epochs. The frequency of the 
observations plotted is 33.6 GHz (black symbols), with the exception of the 2011 June 4 observation, which corresponds to a frequency of 
30.5 GHz (indicated with open symbols). Th e sec ond and third epoch on 2010 November 23 are only separated by ~3 hours. We discuss 
this short-term variability in detai l in Section [321 For source G, the blue open symbols indicate values for extended emission derived from 
smoothed images (see Section [ffTT) . 
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Fig. 5.— Ratio between the flux densities of sources I and BN 
versus time at 33.6 GHz (fllled circles) and at 30.5 GHz (open 
circle). 


gests that the emission from source C is extended, and 
that the higher resolution observations have filtered it 
out. To derive a proper sensitivity level for extended 
emission we have smoothed the B-configuration images 
at Q-band to the resolution of the D-configuration im¬ 
age (1.9"X 1.4"). For the Ka-band monitoring, we have 
smoothed the higher resolution images (from 2011 Feb 
8 to 2011 Jul 09) to a C-configuration resolution of 
0.8"X0.8". We have inspected the smoothed images, and 
measured the flux density (or 3a upper limit) at the loca¬ 
tion of source C. The resulting light curves (Figs. [3] and 
Sj) are in agreement with nearly constant emission during 
the monitoring. This radio source is associated with one 
of the proplyds^^ revealed by the HST (see Fig. E}. Us¬ 
ing the flux densities detected at 33.6 GHz and 43.3 GHz, 
we obtain a spectral index (F ^ z/^) oi a ^3, consistent 
with optically-thin dust emission from the proplyd. 

The other sources detected show clear variations at 
Ka-band between epochs {j3 >0.29). Many of them are 
detected only in some epochs, remaining below the sen¬ 
sitivity limit at other epochs. There is no clear trend or 
pattern in the variability, which appears to be stochastic 
for many of the sources. This highly var iable emission 
sugges ts non-thermal processes. Indeed, iMenten et al.l 
(l2^ detected 4 of these sources (A, F, G and 12) 
with Very Long Baseline Array (VLBA) observations, 
confirming their compactness and hence the non-thermal 
nature of the emission. 

Using epochs separated by ^1 month, it is not possi¬ 
ble to determine whether flux density variation is smooth 
during this period, or whether it happens in shorter 
timescales. Observations with shorter separations are 
needed. We address this issue in Section 13.21 

3.2. Short-term variability 
Day timescales 

The first two Q-band observations are only separated 
by 10 days (2009 March 9 and 19). The main result ob¬ 
tained by comparing these observations is the discovery 
of a new radio source (hereafter OHG-E). Fig. [6] shows 
the comparison of the two images. The source is detected 
at 3.3" northeast of source 1. The flux density variation 
between the two epochs is >5.6 (using the 3cr upper limit 
as the flux density in the 2009 March 9 image). Source 
OHG-E shows a constant flux density during the 4 hours 
of the 2009 March 19 observation, indicating that we 
might have detected a fraction of a powerful flare event. 
This is consistent with the duration reported on other ra¬ 


dio flares observed in Orion (source A and ORBS), which 
have timescales of hours to days. 

Hour timescales 

Two of the Ka-band observations are separated by only 
^3 hours (2010 Nov 23), enabling the study of flux den¬ 
sity variations on even shorter timescales. Eig. [7| shows 
the flux density variation of the 10 sources detected in 
these two epochs, normalized by the flux density of BN. 
In Table [5] we show the ratio between the flux densities 
of the two epochs (F 523 . 89 /F 523 . 76 , where the subscript is 
the JD date minus 2454000, and the variability parame¬ 
ter defined between these two epochs is /3. The sources 
BN, G, 15, 6 are nearly constant within the flux density 
uncertainties, with variations within a factor 0.95—1.17 
and yd < 0.11. The latter 2 sources, however, exhibit 
clear variations at other epochs (see Eig. |4]). 

The sources n, 7, E and 12 show clear radio flares (see 
Eig. [7j), with flux density increases by factors >2.4 and 
variability parameters yd >0.59 (Tab le [5]). The source 
A, a well known radio flaring source (|B^er et al] 120031: 
IZapata et al.l l2004 iGomez et al.l l2008f ). suffers a signif¬ 
icant decrease of its flux density, being undetectable in 
the second image. 

There is a tentative detection of variability in source 
I between the two epochs. The flux densities measured 
are not within the flux density uncertainty limits (see 
Eig. Hand Table[3j), the variability parameter is yd=0.22, 
and the variation factor is 1.37 (while in the case of the 
BN source is only 1.05). We conclude that this short¬ 
term variation of source I is not due to calibration un¬ 
certainties, and might be real. As already commented in 
Sect. 13.11 ioni zed gas of infalling accretion flows onto this 
mass ive star (jGalvan-Madrid et al.lI 2 OIII: iDe Free et al.l 
l2014f ) might be responsible of the long-term variability 
observed in source I (Eig. [5]). However, other mecha- 
nism(s) would be needed to explain the variability ob¬ 
served in hours-timescales, like stoc hastic shocks in th e 
radiative wind of the massive star (jStelzer et al.ll2QQ^ . 
or the presence of an unseen low-mass companion emit- 
ting non-thermal gyrosynchrotron radiation. Indeed, 
iGoddi et al.l (|2Qll[ ) present arguments for source I har¬ 
boring a binary system based on the kinematic history 
of the region. The interaction of a companion with the 
wind from a massive primary could explain the variabil¬ 
ity observed here. 

3.3. Spectral indeces between 30.5 and 37.5 GHz 

One of the methods commonly used to unveil the na¬ 
ture of the radio emission is by studying the emission as 
a function of frequency, F oc z/^, where a is the spectral 
index. 

In the 2011 June 4 observations, we have observed si¬ 
multaneously at two different frequencies, 30.5 and 37.5 
GHz. The calibration uncertainty in this epoch is higher 
than in the other runs (see Section [2]), which implies large 
uncertainties in the derivation of a. However, it is worth 
mentioning that source 7 and source G exhibit clear flux 
decreases with a <-0.4, suggesting that the dominant 
emission mechanism seems to be non-thermal. 


12 p-foplyds are objects for which circumstellar material is being 
ionized by the ultraviolet radiation from massive stars. 


4. COMPARISON WITH LOWER RADIO FREQUENCY 
AND CONNECTION TO X-RAY EMISSION 


























Rivilla et al. 



5''35"'14V4 14V2 


5''35"'14^74 14V2 


RA(J2000) 


Fig. 6.— Upper panels: VLA observations of the OHC region at 45.6 GHz from 2009 March 9 (left) and 2009 March 19 (right). The 
source I is detected in both images, while the new radio flaring source OHC-E appeared 3.3" toward the northeast on 2009 March 19. 
Lower left panel: Zoom-in view of the detection of the source OHC-E from the 2009 March 19 observation. The flux density scale is the 
same as in upper panels. The contours indicate emission from 3(7 to 19(7 in steps of 2(7. Lower right panel: Second detection of OHC-E in 
the 2011 July 9 observation. The contours indicate emission at 3(7, 4(7 and 5(7. 


Table 3 

Multi-epoch radio continuum flux densities®. The values in parenthesis correspond to 3 a upper limits. 


Source 

900 

910 

1188 

1494 

1523.76 

Epoch: JD - 2454000 
1523.89 1569 1601 

1649 

1709 

1717 

1724 

1752 


Q-band (GHz) 





Ka-band (GHz) 






45.6 

45.6 

43.3 

33.6 

33.6 

33.6 

33.6 

33.6 

33.6 

33.6 

37.5 

30.5 

33.6 

33.6 

BN 

31.2T3.3 

30.1±3.1 

27.9±3.0 

26.4T2.7 

25.9±2.7 

27.2T2.8 

23.9T2.5 

24.6±2.5 

25.7±2.6 

27.8±2.8 

20.5i4.1 

18.8i3.8 

25.7i2.6 

22.7±2.3 

I 

14.1±1.7 

12.5±1.5 

13.7±1.7 

9.2±1.1 

9.1±1.2 

12.5±1.5 

10.7±1.3 

10.8il.2 

9.7il.l 

9.3±1.0 

9.1il.9 

7.3il.5 

10.8il.2 

8.1±0.9 

ohc-e 

(1.4) 

7.9±1.1 

(1.6) 

(1.1) 

(1.3) 

(1.2) 

(1.3) 

(1.0) 

(0.7) 

(0.7) 

(0.4) 

(0.3) 

(0.5) 

0.8±0.3 

D 

(1.4) 

(1.1) 

(1.7) 

(1.1) 

(1.2) 

(1.3) 

4.6T0.8 

(1.0) 

(0.7) 

(0.7) 

(0.4) 

(0.3) 

(0.5) 

(0.4) 

n 

(1.4) 

(1.1) 

(1.6) 

4.0T0.8 

(1.2) 

3.5±0.9 

3.5±0.9 

(1.0) 

(0.7) 

(0.7) 

(0.4) 

(0.3) 

(0.5) 

(0.4) 

h 

(1.4) 

(1.2) 

(1.7) 

(1.1) 

(1.2) 

(1.4) 

(1.2) 

(1.0) 

1.4± 0.5 

(0.7) 

(0.4) 

(0.3) 

0.9i0.3 

0.9±0.3 

A 

- 

- 

- 

(5.1) 

9.6±3.8 

(6.2) 

(5.6) 

7.7±3.6 

9.9i2.1 

(3.2) 

- 

9.9i2.4 

(2.2) 

(2.0) 

C 

(35) 

(33) 

27.8T3.7 

12.7T1.7 

11.9±1.9 

11.3T1.9 

12.Oil.8 

19.5±7.2 

20.4±7.4 

15.8±4.0 

(20.8) 

(16.0) 

(18.0) 

(17.4) 

F 

- 

- 

- 

(5.4) 

19.5±4.1 

57.9±7.2 

34.6i5.5 

40.1i5.4 

6.2i2.6 

8.6±2.5 

- 

5.5il.5 

20.0i3.6 

9.9±2.6 

G 

- 

- 

- 

(4.1) 

(4.6) 

(5.2) 

(4.5) 

(3.7) 

(2.5) 

(2.5) 

(2.2) 

5.7il.4 

15.5i2.7 

(1.6) 

e 

(7.1) 

(5.3) 

(6.7) 

(2.3) 

(2.5) 

(2.8) 

7.1±2.4 

11.8±3.1 

(1.4) 

(1.4) 

(1.0) 

(0.5) 

(1.0) 

(0.9) 

15 

(8.5) 

(7.1) 

(8.0) 

(2.6) 

5.9±1.8 

5.6±2.6 

(2.8) 

(2.3) 

3.Oil.3 

(1.6) 

(1.2) 

(0.6) 

(1.1) 

(1.0) 

6 

- 

- 

- 

20.9T3.9 

24.5±4.6 

21.7T4.1 

33.2±4.9 

34.1±5.6 

7.0i2.4 

(3.2) 

- 

16.Oi 3.6 

37.9i6.2 

(2.0) 

7 

- 

- 

- 

(4.1) 

8.0±3.1 

23.8±4.3 

31.4±4.7 

20.7i3.8 

7.1i2.2 

(2.6) 

11.9i3.5 

21.5i4.7 

14.8i 3.6 

(1.6) 

25 

(8.9) 

(7.3) 

(8.2) 

(2.6) 

(2.9) 

(3.1) 

8.8i2.0 

(2.3) 

6.3il.3 

(1.6) 

(1.2) 

(0.6) 

5.4il.2 

19.2±2.1 

12 

- 

- 

- 

(3.1) 

(3.5) 

8.5±2.9 

11.2i2.6 

20.4i2.7 

(1.9) 

38.3±4.2 

18.4i4.3 

13.4i2.7 

69.9 i7.1 

2.8±0.9 

11 

- 

- 

- 

(4.8) 

(5.3) 

(5.9) 

22.8i4.7 

(4.2) 

12.6i3.7 

(2.9) 

- 

27.9i 6.5 

11.9i2.5 

(1.8) 

16 

- 

- 

- 

(5.9) 

(6.5) 

(7.2) 

(6.4) 

(5.2) 

(3.6) 

(3.6) 

- 

(1.1) 

4.7il.6 

(2.3) 

5^ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

30.9i7.5 

- 

- 


“ Integrated flux densities are presented in mJy. 


^ Source 5 is outside the FoV for all but the 2011 Jun 4 at 30.5 GHz, whose FoV is larger. 
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Table 4 

List of radio sources in the ONC/OMC region considered our analysis 


Radio sources 

X-ray sources 


Optical^ 

IR^ 

Mem.^ 

ID 


33.6 GHz 1 

8.3 GHz 


GOUP 


Properties 









Source Z04“ 

F av 

AF 

n 

F av 

AF 

p 

ID 

MedF HRi 

logAn 

logFKs 

BB 

X-ray 

logLx 

OW94 H97 

K05 

HGOO 




(mJy) ( 

mJy) 


(mJy) (mJy) 



(keV) 

(cm“^) 


Num flare?^ 

(erg s-i) 



(propl.) 



I 

19 

10.0 

1.3 

0.13 

0.64 

0.23 

0.36 

/ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

OMG 


14 

14.8 

3.7 

0.25 

4.97 

1.61 

0.32 

- 

- 

- 

- 

- 

- 

- 

- 

y 

y 

y 

ONG 

16 

53 

4.7 

1.7 

0.36 

3.14 

0.43 

0.14 

- 

- 

- 

- 

- 

- 

- 

y 

- 

- 

- 

ONG 

5 

61 

- 

- 

- 

15.88 

0.26 

0.02 

- 

- 

- 

- 

- 

- 

- 

y 

y 

y 

y 

ONG 

BN 

12 

25.5 

1.6 

0.06 

3.61 

0.25 

0.07 

599b'* 

- 

- 

- 

- 

- 

- 

- 

- 

- 

y 

OMG 

OHC-E 

- 

0.8 

0.2 

0.27 

- 

- 

- 

655 

3.58 0.89 

22.85T0.01 

-4.00 

9 

y 

31.37 

- 

- 

- 

y 

OMG 

D 

21 

4.6 

3.0 

0.65 

1.04 

0.68 

0.65 

662 

4.50 0.98 

23.22F0.03 

-4.00 

13 

y 

31.29 

- 

- 

- 

- 

OMG 

n 

17 

3.7 

1.7 

0.45 

1.12 

0.37 

0.33 

621 

3.57 0.84 

22.74F0.01 

-4.00 

13 

y 

30.97 

- 

- 

- 

y 

OMG 

H 

18 

1.1 

0.3 

0.28 

0.66 

0.50 

0.75 

639 

3.99 0.96 

23.06F0.04 

-4.00 

8 

y 

30.90 

- 

- 

- 

- 

OMG 

A 

6 

9.1 

3.7 

0.40 

16.84 

26.7 

1.58 

450 

3.32 0.64 

22.34F0.03 

-4.00 

12 

y 

32.27 

- 

- 

- 

y 

ONG 

F 

76 

24.6 

18.1 

0.74 

15.04 

5.10 

0.34 

965 

1.31 -0.63 

21.25T0.10 

-4.00 

2 

n 

31.89 

- 

y 

- 

y 

ONG 

G 

73 

15.5 

9.8 

0.63 

2.48 

2.09 

0.84 

932 

1.29 -0.64 

21.17F0.09 

-4.00 

12 

y 

32.18 

- 

y 

- 

y 

ONG 

E 

63 

9.5 

5.5 

0.58 

2.07 

0.47 

0.23 

844 

1.63 -0.50 

22.01F0.14 

-0.96 

1 

y 

29.22 

- 

y 

y 

y 

ONG 

15 

49 

4.8 

2.3 

0.48 

4.35 

0.83 

0.19 

766 

1.58 -0.32 

21.54F0.03 

-4.00 

20 

y 

30.71 

- 

y 

- 

y 

ONG 

6 

59 

25.6 

12.8 

0.50 

22.23 

0.33 

0.01 

826 

2.31 0.16 

22.06F0.04 

-4.00 

11 

y 

30.35 

- 

y 

y 

y 

ONG 

7 

52 

17.6 

10.6 

0.60 

10.12 

0.87 

0.09 

787 

2.43 0.29 

22.40F0.04 

-4.00 

5 

y 

30.22 

y 

y 

y 

y 

ONG 

25 

38 

9.9 

6.6 

0.67 

4.96 

0.61 

0.12 

732 

1.33 -0.60 

20.99F1.99 

-4.00 

23 

y 

32.35 

- 

y 

- 

y 

ONG 

12 

41 

25.2 

24.6 

0.98 

12.09 

7.89 

0.65 

745 

1.33 -0.59 

20.79F0.08 

-4.00 

20 

y 

32.33 

y 

y 

- 

y 

ONG 

11 

42 

12.6 

8.5 

0.68 

10.82 

0.17 

0.02 

746 

2.23 0.19 

22.10F0.55 

-3.40 

2 

n 

28.79 

y 

y 

y 

y 

ONG 

ORBS 

- 

- 

- 

- 

- 

- 

- 

647 

5.20 1.00 

23.51F0.03 

-4.00 

8 

y 

30.93 

- 

- 

- 

- 

OMG 

- 

1 

- 

- 

- 

2.84 

0.12 

0.04 

229 

4.05 0.98 

22.96T0.19 

-0.05 

1 

n 

29.28 

- 

- 

- 

- 

EG 

- 

2 

- 

- 

- 

0.81 

0.81 

1.00 

342 

1.98 -0.01 22.21± 0.01 

-4.00 

31 

y 

31.32 

- 

y 

- 

y 

ONG 

- 

7 

- 

- 

- 

0.78 

0.42 

0.54 

510 

4.75 1.00 23.54F 0.06 

-2.59 

3 

y 

31.26 

- 

- 

- 

- 

OMG 

- 

9 

- 

- 

- 

0.48 

0.27 

0.56 

530 

5.23 0.98 

23.50F0.03 

-4.00 

2 

(y) 

30.75 

- 

- 

- 

- 

OMG 

- 

16 

- 

- 

- 

0.33 

0.18 

0.56 

625 

4.81 0.98 

23.42F0.05 

-4.00 

4 

y 

30.95 

- 

- 

- 

- 

OMG 

- 

31 

- 

- 

- 

0.20 

0.08 

0.42 

699 

1.51 -0.50 

21.61F0.il 

-4.00 

4 

y 

29.54 

- 

y 

y 

y 

ONG 

- 

33 

- 

- 

- 

2.90 

1.79 

0.62 

708 

1.54 -0.36 

21.46F0.05 

-4.00 

10 

y 

30.19 

- 

y 

- 

y 

ONG 

- 

34 

- 

- 

- 

5.66 

0.59 

0.10 

717 

1.55 -0.54 

20.98F0.70 

-0.44 

1 

(y) 

28.53 

y 

y 

y 

y 

ONG 

- 

37 

- 

- 

- 

1.44 

0.04 

0.03 

733 

2.49 0.19 

22.54F0.23 

-2.33 

2 

n 

29.40 

y 

- 

y 

y 

ONG 

- 

43 

- 

- 

- 

4.05 

0.43 

0.11 

747 

3.47 0.33 

22.25F0.35 

-1.06 

1 

y 

28.88 

y 

y 

y 

y 

ONG 

- 

44 

- 

- 

- 

1.33 

0.48 

0.36 

757 

1.89 -0.17 

22.45F0.21 

-0.64 

1 

(y) 

29.01 

y 

- 

y 

y 

ONG 

- 

45 

- 

- 

- 

6.03 

0.58 

0.10 

758 

1.66 -0.29 

21.73F0.02 

-4.00 

29 

y 

31.05 

y 

y 

y 

y 

ONG 

- 

46 

- 

- 

- 

1.70 

0.84 

0.49 

768 

1.63 -0.33 

21.52F0.04 

-4.00 

4 

y 

30.27 

y 

y 

y 

y 

ONG 

- 

54 

- 

- 

- 

1.02 

0.61 

0.60 

800 

1.93 -0.14 

22.36F0.21 

-0.85 

2 

y 

28.86 

- 

y 

y 

y 

ONG 

- 

56 

- 

- 

- 

0.65 

0.45 

0.69 

807 

1.37 -0.60 

21.43F0.10 

-4.00 

7 

y 

29.73 

- 

y 

y 

y 

ONG 

- 

58 

- 

- 

- 

1.58 

0.16 

0.10 

820 

3.64 0.95 

22.86F0.13 

-4.00 

3 

y 

30.07 

y 

- 

y 

y 

ONG 

- 

60 

- 

- 

- 

3.12 

0.41 

0.13 

827 

1.63 -0.41 

22.24F0.03 

-4.00 

2 

n 

30.46 

y 

y 

y 

y 

ONG 

- 

64 

- 

- 

- 

7.65 

0.45 

0.06 

847 

1.35 -0.64 

21.62F0.35 

-2.00 

1 

(y) 

29.22 

y 

y 

y 

y 

ONG 

- 

65 

- 

- 

- 

3.82 

0.18 

0.05 

855 

1.96 -0.01 

21.69F0.09 

-4.00 

8 

y 

31.79 

y 

y 

y 

y 

ONG 

- 

69 

- 

- 

- 

2.03 

0.38 

0.19 

876 

3.22 0.70 

22.69F0.06 

-2.60 

3 

y 

29.81 

- 

y 

y 

y 

ONG 

- 

71 

- 

- 

- 

4.05 

0.39 

0.10 

900 

4.04 0.46 

23.54F0.46 

-4.00 

2 

y 

30.44 

- 

y 

y 

y 

ONG 

- 

75 

- 

- 

- 

0.43 

0.14 

0.33 

955 

1.70 -0.31 

21.46F0.19 

-0.36 

2 

y 

28.78 

- 

y 

y 

y 

ONG 

- 

77 

- 

- 

- 

2.81 

1.71 

0.61 

1130 

1.23 -0.71 

21.28F0.14 

-4.00 

7 

(y) 

31.67 

- 

y 

- 

y 

ONG 

- 

3 

- 

- 

- 

0.16 

- 

- 

378 

1.23 -0.66 

21.11F0.07 

-4.00 

5 

y 

30.20 

- 

y 

- 

y 

ONG 

- 

4 

- 

- 

- 

0.30 

- 

- 

394 

1.20 -0.72 

20.00F1.65 

-4.00 

6 

y 

31.32 

- 

y 

- 

y 

ONG 

- 

5 

- 

- 

- 
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^ Source ID from f^apata et al.l 1200411 . 

^ Optical association from the work of H97 IHillenbrandll 19971 ') and K05 IKastner et aT1l200^ '). 

^ Infrared association from the work of HCOO m illenbrand &: Carpenterll200^~) . 

^ Source membership: OMC = Orion Molecular Cloud, ONC = Orion Nebula Cluster, EG = extragalactic. 

® From visual inspection of the X-ray light curves published bv IGetman et al.l 120051') . The parentheses denote that the detection is uncertain. 


Although X-ray emission from wind shocks would be expected from the massive star associated with source I, the non-detection by Chandra is likely 
presence of a nearly edge-on disk IMatthews et al.1120101 ') that absorbs the X-ray emission. 


As discussed in Section 13.11 the source C is filtered out by the more extended VLA configuration observations. To compute the average flux, the 
parameter and the number of detections at 33.6 GHz, we have considered the images after smoothing to a C-configuration resolution of 0.8^^. 


due to the 
variability 


^ The massive BN object has an X-ray counterpart, GOUP 599b, much fainter than a low-mass companion located 0.9^^ from BN, COUP 599a IGrosso et al1l2005l '). 
h 


The radio source 13 identified bv IZapata et al.l 120041 ) does not have optical or IR counterpart, but it has an X-ray counterpart (COUP 607) with low 
counts, which prevent the derivation of Npr. IRivilla et al.l 120131) showed that this sources is likely driving a molecular outflow in the OMCl-S region. 


number of 
and hence 
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Julian Date — 2454000 

Fig. 7. — Flux density variation (normalized by the flux density 
of BN) on hour timescales for sources detected in the two runs 
observed on 2010 November 23. 



30 35 


V (GHz) 

Fig. 8.— Flux densities at 30.5 and 37.5 GHz from the 2011 
June 4 observations. 


With the aim of better understanding the nature of 
the radio emission from young stars of the ONC/OMC 
region, in this section we compare the properties of the 
33.6 GHz emission from o ur monitoring with those of 
the 8.3 GHz emission from IZapata et aO (|2QQ4l ). These 
authors analyzed VLA observations in 4 different epochs 
within a larger FoV and cadence year. We have cross- 
correlated our sample of 18 radio sources emitting at 33.6 
GHz with their catalog of radio sources. 

Furthermore, to explore the link between the radio and 
X-ray emission in a statistically significant sample, we 


Table 5 

Radio variability on hour timescales between the 
TWO OBSERVATIONS ON 2010 NOV 23. 


Source 

^523.89/^523.76 


BN 

1.05 

0.03 

I 

1.4 

0.22 

n 

2.9 

0.69 

F 

3.0 

0.70 

7 

3.0 

0.70 

A 

0.7 

0.30 

G 

0.95 

0.04 

15 

0.95 

0.04 

6 

1.17 

0.11 

12 

2.4 

0.59 



Fig. 9. — Scheme of the method used to classify the sources as 
ONC or OMC members, or EG sources. We use cross-correlation 
with optical, IR and X-ray stellar catalogs, and comparison with 
the spatial distribution of molecular material (see text). 


study the full sample of sources emitting at both wave¬ 
lengths in Orion. We have cross-correlated the full sam¬ 
ple radio sources with the catalog of X-ray stars provided 
by the very deep Ghandra Orion Ultra Deep Project 
(COUP, iGetman et al.ll2QQ5l) . 

The ONC/OMC region harbors different populations 
of radio stars: i) members of the optically visible fore¬ 
ground ONC illuminated by the massive stars of the 
Trapezium (some of them associated with proplyds); ii) 
stars still embedded in the background OMC in a earlier 
evolutionary stages; and iii) extragalactic (EG) sources. 
Since the emission mechanisms of each group could be 
different, we distinguish in our study these different pop¬ 
ulations. 


4.1. Source memhership 

To obtain a rough e stimate of the number of expected 
EG sources we follow Ihbrnalont et al.l (|2QQ2[) . who esti¬ 
mated the expected number of EG contaminants at 8 
GHz. In a typical EoV at 8 GHz, the expec ted EG 
contamination is source ([Zapata et al.ll2QQ^ . In the 
smaller EoV of our 33.6 GHz monitoring, we would ex¬ 
pect an even lower EG contamination of ^0.01 sources at 
8 GHz for the most sensitive of our observations and a de- 
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Fig. 10.— Spatial distribution of the full sample of sources used 
in our analysis. Red and blue dots indicate members of the ONC 
and OMC, respectively. The ONC sources related with proplyds 
are denoted by large blue circles. The radio source classified as EG 
falls outside the region shown. The greyscale-scale shows the IR 
K-band image from the 2 Micron All Sky Survey (2MASS). The 
dashed contours trace gas from the OMC (CN N=l—0 emission 
from [Rodriguez-Franco et al.lfTOOSIb from 12 K km s‘^ in steps of 
4 K km s'^. The dotted contours correspond to emission at 850 
p,m from IDi Francesco et al.l (|2008l b The large circle indicates the 
primary beam of the 33.6 GHz observations. 


tection criterion of 5cr. Since these extragalactic sources 
generally show non-thermal emissio n of the type F oc 
with (a < 0 between 1 — 100 GHz (| Condon! 1 1992f ). the 
contamination at 33.6 GHz is expected to be even lower. 
We therefore conclude that the impact of EG contam- 
inat ion in onr m o nitori ng is very low. In the case of 
the IZapata et al.l (|2004f ) observations at 8.3 GHz, they 
inferred ^1 EG source. 

We classify the radio sources into 4 groups: ONC mem¬ 
bers without evidence of proplyds (“naked” ONC), ONC 
members associated with proplyds, embedded OMC 
members and EG contaminants. We cross-check the ra¬ 
dio sources with available stellar catalogs in the opti- 


__ _ log s _ . 

cal (|0’dell fc Wenlll994lHillenbrandlll9971:lKastner et al.l 

|2QQ5[) . infrar ed (iHillenbrand fc Caroenterl l2QQQ[) and X- 

rays (COUP. FGetman et al.ll2QQ5f ). Additionally, we com¬ 
pare the position of the sources with the spatial distri- 
butio n of the OMC, traced b y the C N N=l—0 emission 
from iRodriguez-Eranco et al.l (|l998f ). A scheme explain¬ 
ing the method used to classify the sources is presented 
in Eig. [9l It is based on 6 steps: 


1. The presence of optical counterpart indicates that 
the source is likely an ONC member. Since it is 
possible that some of the sources identified in this 
way are foreground field stars not related with the 
young cluster, we have additionally cross-correlated 
the radio sample wit h the list of 16 field stars from 
iGetman et al.l (|2QQ5f ). without any coincidence. 

2 . The presence of IR counterpart indicates that the 


source is likely an ONC or OMC star, because EG 
sources are expected to be weaker IR sources. 

3. The value of hydrogen column density A^h derived 
from X-rays is in general a good indicator to dis¬ 
criminate between ONC and OMC members. A 
source with logA^n <22.5 cm“^ is considered as 
ONC member 

4. The presence of X-ray variability is a good indica¬ 
tor to determine if a highly extincted source with¬ 
out IR counterpart is an embedded OMC member 
or an EG source, because young stars are expected 
to exhibit much h i gher X -ray variability. Eollow- 
ing IGetman et al.l (|2QQ5t ) we consider 3 signposts 
of X-ray radio variability: i) the significance of a 
Kolmogorov-Smirnov test (Pks), which establishes 
whether variations are above those expected from 
Poisson noise associated with a constant source; 
ii) the number of segments of the Bayesian block 
parametric mod el (BBNum) oi source variabil¬ 
ity developed by IScarglel (|l998f ) ; hi) visual inspec¬ 
tion of the X-ray light curves. We consider that a 
source is X-ray variable when Pks < — 2.0 and/or 
BBNum > 2 and/or it exhibits X-ray flares in 
the light curves. We have obtained the values of 
Pks and BBNum, and visually examine d the X- 
ray light curves from IGetman et ahl (|2QQ5t ) (see Ta¬ 
ble |4]). 

5. Einally, we compare the position of the sources with 
respect to the location of the OMC, traced by the 
CN N =l—0 emission from iRodriguez-Eranco et al.l 

(see contours in Eig. [iQl). 

Additionally, to identify the ONC stars related to pro¬ 
plyds, we have cross-correl ated the radio samp le with 
the catalog of proplyds from iKastner et aP (|2QQ5f ). using 
a counterpart radius of 0.5". The resulting membership 
classification is shown in the last column of Table IH Eig. 
M shows the spatial distribution of the different groups 
of sources. 

4.2. Sample of sources emitting at 8.3 and 33.6 GHz 

The cross-correlation between our radio sample at 33.6 
GHz and the radio sample at 8.3 GHz, with a counterpart 
radius of 0.5", shows that all the 33.6 GHz sources were 
also detected at 8.3 GHz, with the only exception of the 
new radio flaring source OHC-E. On the other hand, our 
monitoring did not detect emission from 25 sources from 
IZapata et al.l (|2004l ) located within our EoV. The radio 
properties at both radio frequencies are summarized in 
Table IH In the left panel of Eig. [TT] we compare the 
average flux densities at both frequencies. We have dis¬ 
tinguished “naked” ONC stars, ONC stars related with 
proplyds and OMC stars. Eor those sources without 33.6 
GHz detection, we have considered 3cr upper limits mea¬ 
sured in our image with best sensitivity. 

In general, the radio sources detected by our moni¬ 
toring exhibit higher 8.3 GHz flux densities than those 
undetected. This suggests that the non-detection at high 
frequency is likely due to lack of sensitivity. 

With only one exception, the radio sources have higher 
average flux densities at 33.6 GHz. In the case of sources 
BN, I and C, we interpret that this is due to “quiescent” 
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Fig. 11.— Comparison of radio properties at 33.6 and 8.3 GHz. Red and blue dots indicate members of the ONC and CMC, respectively. 
The ONC sources related with proplyds are denoted by large blue circles. Left panel: Mean flux density at 33.6 GHz versus the mean flux 
density at 8.3 GHz. For those sources that fall within the FoV of our 33.6 GHz observations but are only detected at 8.3 GHz, we have 
included 3f7 upper limits (blue solid triangles; and large empty triangles if they are related with proplyds). The dashed lines show the 
relation expected for different types of emission: optically-thin thermal dust (spectral index a ~ 3); optically-thick, ionized, stellar wind 
{a = 0.6); optically-thin free-free (a = —0.1); and non-thermal synchrotron (a = —0.7). We also indicate with a solid line the case of a 
ffat spectrum (a = 0). Middle panel: Comparison between AF at 8.3 GHz and 33.6 GHz. The symbols are the same as in the left panel. 
Right panel: Comparison between the variability parameters (3 at 8.3 GHz and 33.6 GHz. The symbols are the same as in the left panel. 


8.3 GHz 


33.6 GHz 




6av (mJy) (mJy) 

Fig. 12.— Radio properties at 8.3 GHz (left) and 33.6 GHz (right) of the radio/X-ray sample. The black diamond corresponds to the 
EG source. The other colors and symbols are the same as in Fig. (TO] The dashed lines indicate values of /3 of 1, 0.1 and 0.01. 


thermal component that increases with frequency. In the 
other sources, this does not directly imply evidence of 
thermal origin, because the observations are not simul¬ 
taneous and their emission is highly variable (Section [3j). 
Since this variability is likely connected to non-thermal 
processes, the value of Fav is very sensitive to the pres¬ 
ence of flares, and not representative of their thermal 
emission. 

In the middle and right panels of Fig. [TT]we compare 
the values of AF and (3 at 8.3 and 33.6 GHz. Most of 
the sources appear more variable at higher frequencies 
(15/17 have higher AF; and 11/17 have higher /3). This 
may indicate that the radio variability increase with fre¬ 
quency, although new observations of a more statistically 


representative sample of sources are needed to confirm 
this behavior. 

We note that even the sources related with proplyds, 
which are expected to emit non-variable free-free or dust 
emission from circumstellar material, are clearly variable. 
This points toward the presence of a non-thermal com¬ 
ponent arising from the central PMS stars. 

4.3. Full sample of sources emitting radio and X-rays 

We have cross-correlated of the sample radio sources 
(detected at 8.3 GHz and/or at 33.6 GHz^^) with the 
COUP catalog, searched for X-ray counterparts within 

We have also includ ed the ffaring radio source ORBS detected 
bv IForbrich et aP (|2008l f at 22 GHz. 
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0.5" of the radio sources. We have obtained a final sam¬ 
ple of 51 sources emitting at radio and X-rays. The prop¬ 
erties of the X-ray counterparts are shown in Table IH 

Radio properties 

In Fig. [12] we show the radio properties of the subsam¬ 
ples of X-ray sources with emission at 8.3 and 33.6 GHz. 
The ONC sources related with proplyds exhibit low ab¬ 
solute variability at 8.3 GHz {AF < 1 mJy) compared to 
“naked” ONC stars, which appear more variable (AF > 
1 mJy). This suggests that the 8.3 GHz emission in the 
proplyds is dominated by thermal nearly constant emis¬ 
sion from circumstellar material, while “naked” ONC 
stars are might be dominated by variable non-thermal 
processes related with magnetic activity of the PMS star. 
This could be gyrosynchrotron emission produced by the 
acceleration of electrons in mag netic field re connection 
events in the corona of the star ([Andrei 1 19960 . 

However, we note that most of the proplyds with lower 
Fav have values of P >0.1 (left panel of Fig. [T^ . because 
their AF is a significant fraction of their Fav There¬ 
fore, it seems that the 8.3 GHz emission from proplyds, 
although dominated by thermal emission, could include 
also a non-thermal component arising from the central 
PMS stars. 

The OMC members show low values of absolute vari¬ 
ability at 8.3 GHz (AF < 1 mJy). However, since most 
of them are weak sources, this low absolute variation 
translates into significant relative variability, with values 
of yd >0.1, also pointing to non-thermal origin. Only BN 
is nearly constant, with low values of both AF and yd, as 
occurs at 33.6 GHz. This confirms that the radio emis¬ 
sion is likely dominated by thermal emission from the 
ionized gas in the massive envelope around the star. 

The right panel of Fig. 0 shows that all sources appear 
in general more variable at 33.6 GHz (both in AF and yd). 
This may indicate that the emission at higher frequen¬ 
cies is more dominated by non-thermal highly variable 
processes. 

Therefore, we conclude that the radio emission can 
be attributed to two different mechanisms: i) highly- 
variable (flaring) non-thermal radio gyrosynchrotron 
emission produced by accelerated electrons in the magne¬ 
tospheres of low-mass PMS stellar members of both the 
ONG and OMC; and ii) non-variable thermal emission 
from the ionized gas and heated dust of the ONC pro¬ 
plyds illuminated by the Trapezium Cluster, or from the 
ionized gas in the envelopes surrounding massive stars, 
as in the case of the BN object. 

X-ray properties 

• Hydrogen eolumn density Nn, Hardness ratio (HRi) 
and Median energy (MedE): in Fig. [13] the values of 
median energy {MedE) of X-ray photons and the hard¬ 
ness ratio^^ {HRi) of the full radio/X-ray sample as a 
function of Ah are shown. It is clear that the stars em¬ 
bedded in the OMC have higher v alues (MedE >3 . 5 keV 
and HRi >0.84) than ONC stars. iFeigelson et al.l (|2QQ5[ ) 
found a relation between the median energy of the X-ray 

Hardness ratio HRi = (h — s)/(hs)^ where h and s refer to 
the counts detected in the hard (2.0—8.0 keV) and soft (0.5—2.0 
keV) bands, respectively. Values closer to 1.0 indicate a hard X-ray 
source, and —1.0 a soft X-ray source. 


photons arising from the stars and the hydrogen column 
density Ah (left panel of Fig. [T3|) . This trend is due to 
an absorption effect: in the embedded sources, only the 
harder photons are able to escape through the molecular 
gas and then be detectable, while softer ones (with lower 
energy) are absorbed. As a consequence, the embedded 
stars appear as harder sources. Among the ONC mem¬ 
bers, those related with proplyds exhibit higher values of 
Ah (and hence of MedE and HRi). This is likely due 
to the presence of circumstellar material that produces 
higher extinction than in “naked” ONC stars. 

• X-ray luminosities: In Fig. [IT] we show the ra¬ 
dio properties of the radio/X-ray sample as a function 
of the extinction-corrected X-ray luminosity Fx- We 

find that lONC-proplyds ^ j^oMC ^ j^ONC-naked^ 

though this result could be physically real, we note that 
it should be taken with caution, because the derivation 
of the extinction-corrected luminosity is more uncertain 
for highly extincted sources. 

We also investigate whether there is a relation between 
the radio properties and the X-ray l u minos ity in the 
radio/X-ray sample. iForbrich fc Wollj (|2Q13[) remarked 
that there is no clear c orrelation between t he radio flux 
densities measured by IZapata et al.l (|2QQ4l ) at 8.3 GHz 
and the COUP X-ray luminosities (upper left panel in 
Fig. [T4|). Similarly, we do not see a trend in our subsam¬ 
ple at 33.6 GHz (upper right panel in Fig. [T4|) . How¬ 
ever, if we do not consider the proplyds a tentative trend 
appears: the radio flux density generally increases with 
the X-ray luminosity. This suggests a relation between 
the X-ray and radio emissions for the “naked” ONC and 
OMC sources. In the case of proplyds, a significant frac¬ 
tion of radio emission is expected to arise from gas ion¬ 
ized by external illumination, which is not linked to X-ray 
activity. Then, a relation between the radio flux and X- 
ray luminosity is not expected in proplyds, as observed 
in Fig. [14] 

We have also studied the fraction of X-ray sources de¬ 
tected at radio wavelengths as a function of the X-ray 
luminosity (Fig. [T5|) for the 8.3 and 33.6 GHz monitor¬ 
ings. In each case, we have considered the COUP sources 
that fall within the FoV of each observation with enough 
counts so that the X-ray luminosity corrected for extinc¬ 
tion could be reported. We find 159 X-ray sources within 
our FoV and 595 sources within the 8.3 GHz FoV. Both 
samples show a very similar behavior. The weaker X-ray 
sources with logFx < 28.5 erg s“^ (200 sources, i.e., 13% 
of the full COUP s ample) were n o t dete cted by our radio 
monitoring, while IZapata et al.l (|2004l ) only detected 1 
source. Fig. [15] shows that the fraction of X-ray sources 
detected in the radio increases with X-ray luminosity. 
This indicates that the radio observations have statisti¬ 
cally detected the most luminous X-ray sources. This 
suggests that: i) the underlying mechanisms responsible 
for the X-ray and (at least some fraction of) the radio 
emission are somehow related; and ii) radio monitorings 
have been limited in the past due to sensitivi ty, confirm¬ 
ing the conclusion of iForbrich X WoTfl (|2Q13D . 

The middle and lower panels of Fig. [14] also show ten¬ 
tative trends between the radio variability (AF and /3) 
and the X-ray luminosity for the sources not related with 
proplyds. In general, the radio variability increases to- 
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Fig. 13. — Relation between the hydrogen column density Nu and the median energy of X-ray photons {MedE, left panel) and the 
hardness ratio ( if Ri , right panel) of the members of the radio/X-ray sample. The dashed line in left panel is the empirical fit found by 
IFeigelson et al.l (I2005I V The colors and symbols are the same as in Fig. 1121 


wards higher Lx- 

However, we note that the number of sources of the 
“naked” ONC and OMC subsamples is too low to draw 
robust general conclusions. Deeper radio observations 
detecting a much larger number of sources are needed to 
confirm these tentative trends. 

• X-ray variability: We now study whether there is 
a relation between the radio and the X-ray variability. 
We follow the 3 conditions presented in Sect. 14.11 to 
consider a source X-ray variable: Pks < —2.0 and/or 
BBNum > 2 and/or presence of X-ray flares in the 
light curves (see Table |4]). Only one radio source (E) of 
the 33.6 GHz/X-ray subsample does not exhibit X-ray 
variability. Namely, 93% of the stars of this subsample 
are X-ray variable. Regarding the full radio/X-ray sam¬ 
ple, 85% of the stars are X-ray variable. These fractions 
are higher than the fraction of X-ray vari able sources of 
the f ull COUP sample, which is ^60% (|Getman et al.l 
l2QQ5f ). Therefore, we have found that the radio sources 
are mostly associated with X-ray variable stars. This 
supports that (at least some of the) radio emission might 
be related to the same magnetic events in the coronae of 
PMS low-mass stars that also produce the X-ray emis¬ 
sion. 

In Fig. [16] we show the radio properties of the radio/X- 
ray sample as a function of the level of X-ray variabil¬ 
ity, quantified by BBNum. As in Fig. [HI tentative 
trends appear if we do not consider the proplyds. How¬ 
ever, given that the X-ray and radio observations were 
carried out at different epochs and different sensitivities, 
no firm conclusions can be drawn. Unlike the very deep 
COUP observation, which continuously observed the re¬ 
gion during ^10 days, the radio monitoring runs are 
much shorter, yielding a lower probability to detect large 
flux density variations such as flares. To better under¬ 


stand if radio and X-ray variability are directly related, 
simultaneous observations will be needed. 

5. THE NEW RADIO SOURCE EMBEDDED IN THE ORION 
HOT CORE: OHC-E 

Our observations have detected twice a new radio 
source, OH C-E (2009 March 1 9 and 2011 July 09, s ee Fig . 
[6|). Neither iFelli et al.] (|l993f ) nor IZapata et al.l (|2004 ). 
who performed observations covering 7 months and 4 
years at 5 and 15 GHz and 8.3 GHz, respectively, nor 
iGoddi et al.l (|201ll ). who observed the same region at 45 
GHz two months before our detection (2009 January 12), 
detected emission towards this source. 

The source is not resolved by the beam of our observa¬ 
tions in any of our two detections (^ 0.2" and ^ 0.075"). 
We can set upper limits to the deconvolved size of the 
emission of < 0.1" (< 40 AU) for the Q-band emission 
and of < 0.04" (< 10 AU) for the Ka-band emission. Un¬ 
fortunately, it is not possible to determine the spectral 
index of the emission, because observations at several 
radio wavelengths were not carried out simultaneously. 
We derived lower limits to the brightness temperature 
of ^ 460 K and ^ 615 K from the source flux density 
at 45.6 GHz and 33.6 GHz, respectively, and the upper 
limits to the source sizes. Although these temperatures 
are consistent with both thermal and non-thermal emis¬ 
sion, our monitoring shows that the emission is highly 
variable, suggesting a non-thermal origin. 

The position of the radio flare coincides with an em¬ 
bedded X-ray low-mass pre-main sequence (PMS) star 
COUP 655 (Fig. [Hj). The source OHC-E is also located 
very close to the southeast member of a binary stellar 
system (CB4, Fig. [T7|) observed with high angular res¬ 
olution by the Near Infrared Camera and Multi-Object 
Spectrometer (NICM OS) onboard the Hubble Space tele¬ 
scope (HST) (see also lStolovv et all! 1 9981: ISimpson et al.l 
|2QQ6[ ) . Therefore it seems that both the radio and the 
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Fig. 14.— Radio properties (Fav, AF, /3) at 8.3 and 33.6 GHz versus the value of the X-ray luminosity derived from the X-ray COUP 
counterparts. The colors and symbols are the same as in Fig. 1121 
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Fig. 15.— Fraction of X-ray COUP sou rces detected by our 
radio monitoring at 33.6 GHz (solid line) and lZaoata et all (120041) 
monitoring at 8.3 GHz (dashed line) as a function of the absorption- 
corrected total X-ray luminosity Lx- In each case, we have consid¬ 
ered the COUP sources that fall within the FoV and with enough 
counts so that the X-ray luminosity corrected for extinction could 
be reported: 15 9 for our monitoring and 595 for the 8.3 GHz mon¬ 
itoring (see also IForbrich fc Wolk|[20T^) . 


X-ray emission are related to this star. We can give a 
rough estimate for the mass of this star using the canon- 
i cal relation for X-ray young stars, log[I/x/ivboi] =—3.0 
(|Pallavicini et al.l[l98lf ). Assuming an stellar age typical 
for a massive star forming region of 5 x10^ yr, OHC-E 
would have a mass of Mq (using the lSiess et al iiinno 
stellar models), supporting that it is a low-mass star. 
Hence, the variability observed both in radio and X-rays 
wavelengths can be related with magnetic activity in the 
corona of this PMS low-mass star. 

6. SOURCE 12: THE BINARY 0^ ORI A 

0^ Ori A is a well known binary system, with a BO.5 
prim ary and a low-mass companion (e. g.. IClose et ah 
2013 ). As observed in the binary WR140 (jwilliams et al. 
1990[ ). the binarity can produce a smooth periodic vari¬ 
ation of the radio emission caused by the combination of 
two geometrical effects: i) variation of the free-free opac¬ 
ity due to the ionized envelope of the primary star as the 
companion orbits; ii) variation of the stellar activity (and 
hence the non-thermal emission) inversely proportional 
t o the separation between components. 

iFelli et ahl (|l99lL 1 19931 ) monitored the radio emission 
from this binary at 5 and 15 GHz. We plot in Fig. [18] 
their flux densities and those of our monitoring as a func¬ 
tion of the orbital phase 0. We have used th e orbital pa¬ 
rame ters P=65.4325 and To=JD 2446811.95 (|Bossi et al.l 
Il989[ ). The flux density at 5 and 15 GHz peaks near peri- 
astron (0 ^0.15), with lower levels at (j) ^0.1 and (j) =0.6 
to 0.9. Our higher frequency flux densities are consistent 
with those at lower fr equencies, with a peak after peri- 
astron at 0 ^0.2. As iFelli et aP (|l993f ) noted and our 
data at higher frequency confirm, the orbital modula¬ 
tion model may make the main peak be detected always 
at the same position 0 ^0.15—0.2. However, the pres¬ 


ence of large scatter in t he radio ein i ssion cannot be ex¬ 
plained with this model. iFelli et ahl (| 19931) set an upper 
limit for the variability timescale of 10 to 20 days. Our 
monitoring has revealed variability at shorter timescales 
of hours (Section l3.2[) . This is in agreement with non- 
thermal emission due to stellar activity, perhaps in addi¬ 
tion to orbital modulation. Therefore, we conclude that 
although certain orbital modulation may be present, pro¬ 
ducing the observed flux density peak, there is also a 
non-thermal emission component arising likely from the 
low-mass companion that varies independently of the or¬ 
bital phase. 

7. RATE OF RADIO FLARING AGTIVITY 

There are only a few radio flares detected from young 
stellar objects. This could indicate that these phenom¬ 
ena are rare events, or alternatively, that the sensitivity 
of the observations carried out so far has limited their 
detection. Our monitoring allows us to evaluate which 
is the most likely scenario. We have detected flares in 
sources OHG-E, F, n, 7 and 12. Other sources (e.g. D, 
E, G, 6, 11, 12 and 25) showing high levels of variability 
through different epochs separated by longer timescales, 
might be also flaring sources, but with our data we are 
n ot able to con f irm sh ort-term variability. 

iBower et al.l (j2QQ3f ) parametrized the rate of radio 
flares (flares day“^) in the entire ONG/OMG region 
based on their d ata and the previo us observations of vari¬ 
able sources by lEelli et ^ (|l993f ) . Since the number of 
flares depends on the field of view (EoV) of a particular 
observation, and of the stellar density of the observed re¬ 
gion, we reformulate the expression for the rate of radio 
flares, A^rf, to: 


A^rf = 7 


100 mJy 


a 


f ^FoV 

\ ^tot 



( 1 ) 


where F^a is the threshold detection limit (considered 
as 3 times the RMS of the observation); a is the spec¬ 
tral index of the emission; ApoV and Atot are the areas 
of the observed region and the full ONG/OMG region, 
respectively; and EpoV and E are the surface stellar den¬ 
sity in the observed region and the mean stellar density 
of the entire ONG / OMG, respectively. The parameter 
7 is a constant that iBower et al.l (|2003f ) estimated to be 
between 0.01 and 0.1. 

Since the flare rate is a function of the stellar density, 
it is convenient to observe the more crowded region of a 
cluster to enhance the chances of detecting radio flares. 
The r egion with highest stellar density in Orion is the 
OHG (|Rivilla et al.ll20l'^ . Gonsidering a ^25"x25" re¬ 
gion centered in the stellar density peak, using the typical 
power-law index of non-thermal emission of a = — 1, the 
mean RMS noise of our monitoring (0.3 mJy), the ap¬ 
proximated size of the ONG/OMG cluster of 15 pcxl5 
pc, and the census of stars in the region provided by 
GOUP, we obtained that A^rf ^0.017-0.17 flare day“^. 

We have detected 2 clear flares (OHG-E and source 
n^^.) in this area in 14 observations, so a rough es¬ 
timate of the flaring rate is ^0.14 flare observation”^. 


The detections of sources H and D in the OHG region show 
also ten tativ e evidences of flaring emission towards these stars (see 
Section Ha 
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Fig. 16.— Radio properties (Fav, AF, /3) at 8.3 and 33.6 GHz versus the value of the BB Num derived from the X-ray COUP counterparts, 
which is a proxy for the X-ray variability. The colors and symbols are the same as in Fig. 1121 
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Fig. 17.— 2.15 /j,m HST-NICMOS image (color scale) of the re¬ 
gion where OHC-E was detected, from the Hubble Legacy Archive 
(HLA). An IR binary is detected toward this region, CB4 (see also 
IStolovv et al. I [19^ : I Simpson et al.ll2006lL The green contours cor¬ 
respond to the radio emission at 45.6 GHz detected on our 2009 
March 19 image (5cr, lOcr and Ihcr). The positions of the X-ray 
star COUP 655 is indicated with the white plus sign. 


GHz (band 3) with a full bandwidth of 7.5 GHz and 50 
antennas, can reach a 8 /iJy sensitivity limit in only 2.3 
hours of on-source observing time^^ Gonsidering a FoV 
of ^25"x25", this ALMA observation may find ^ 6 ra¬ 
dio flares day“^, which represents ^ 25% of the X-ray 
sources detected by Ghandra in the region. This would 
confirm that radio flares are common events, similarly to 
the X-ray flares detected by Ghandra. 

Future observations are clearly needed to derive a bet¬ 
ter estimate of this radio flaring rate, through the detec¬ 
tion of many more radio flares. The flaring information 
for tens or even hundreds of PMS stars will provide a 
complete statistical description of radio short-term vari¬ 
ability. Beyond the purely scientific interest, this would 
have imp ortant technical irn plications for interferometric 
imaging (|Bower et all 120031) . The classical interferomet¬ 
ric imaging techniques assume a constant sky in image 
reconstruction. However, this assumption would be vi¬ 
olated by the presence of many variable sources in the 
field. This would lead to a reduced dynamic range of the 
image ([Stewart et al.ll2QTlh . Also, it would become dif¬ 
ficult to concatenate multiple observations of the same 
region to obtain deeper images. Therefore, a more accu¬ 
rate determination of the radio flare rate would help to 
understand to what extent this may affect deeper VLA 
and ALMA observations of young stellar clusters. 



Fig. 18.— Flux densities of source 12 as a function of the orbital 
period at 33.6 GHz (red dots, this paper), 5 a nd 15 GHz (blue 
squar es and green triangles, respectively, from IFelli et ahl 119911 . 

\im)- 


Assuming that the typi cal duration of the radio flares 
is several hours to d ays ([Andrei[T q 9 61: iBower et al.ll2QQ3l: 
iForbrich et al.ll2QQ^ . this would be approximately equiv¬ 
alent to 0.14 flare day“^, which is similar to the upper 
value from Eq. [TJ Therefore, our results suggest that 7 is 
closer to 0.1 rather than to 0 . 01 , and consequently that 
the number of detectable flares is significant. Therefore, 
our multi-epoch monitoring confirms that the presence of 
radio flares is not a rare phenomenon in crowded young 
stellar clusters. 

Obviously, the detection of radio flares in the OHG is 
favored by its high density of embedded low-mass stars. 
But, according to Eq. [H even in less dense regions the 
number of detectable flares would be significant, espe¬ 
cially if the sensitivity is enhanced. The improved capa¬ 
bilities of the VLA and ALMA may be expected to reveal 
many more radio flares arising from young low-mass star 
clusters. A single polarization ALMA observation at 90 


8. SUMMARY AND GONGLUSIONS 

In this work we have presented a multi-epoch radio 
monitoring of the ONG/OMG region carried out with 
the VLA at high centimeter frequencies (33 and 45 GHz). 
We have detected 19 radio sources, mainly concentrated 
in the Orion Hot Gore and Trapezium regions. Two of 
them are related with massive stars: sources BN and 1. 
The flux densities of the BN source and source G (re¬ 
lated with a proplyd) are compatible with constant ther¬ 
mal emission. The source I, besides a constant thermal 
component, shows tentative evidence of radio variability 
at both short- and long-term timescales. The remaining 
16 sources show long-term (month-timescale) variability, 
but it is not yet clear whether these comprise multiple 
short-term events not covered by our monitoring cadence. 
Indeed, we have confirmed radio flares (i.e., short-term 
radio variability on timescales of hours to days) in 5 
sources: E, 7, n, 12 and the new source OHG-E, pre¬ 
viously undetected at radio wavelengths. 

We have complemented our radio sample with other ra¬ 
dio detections at 8.3 GHz from the literature, and cross- 
correlated it with the X-ray GOUP catalog to obtain the 
full sample of sources emitting radio and X-ray in the 
ONG/OMG region. The radio emission from young stars 
can be explained by a combination of 2 different mech¬ 
anisms: i) non-variable thermal emission produced by 
ionized gas and/or heated dust from the ONG proplyds 
and the massive objects BN and I; and ii) variable (flar¬ 
ing) non-thermal gyrosynchrotron emission produced by 
accelerated electrons in the stellar corona of PMS low- 
mass members of the ONG and OMG. We have found 
several hints relating this variable radio emission with 
the X-ray activity. 

Our study of the radio variability of 0^ OriA concludes 

According to the ALMA sensitivity calculator available in the 
ALMA Observing Tool. 
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that there is evidence of a non-thermal emission compo¬ 
nent arising likely from the low-mass companion. More¬ 
over, certain orbital modulation may be present in this 
binary, producing the observed flux density peak. 

We have derived a rough estimate of the radio flaring 
rate in the densest cluster in the region, which is em¬ 
bedded in the Orion Hot Core. We have obtained ^0.14 
flares day“^. This value is consistent with a empirical es¬ 
timate assuming the sensitivity and FoV of our observa¬ 
tions and the stellar density of the region. This confirms 
that radio flares are not rare phenomena during the ear¬ 
liest stages of star formation as previously thought, but 
relatively common events similarly to the well-known X- 
rays flares. 

Our results have shown that the radio monitorings to 
date has been strongly limited by sensitivity, detecting 
mainly those sources with higher X-ray luminosity. This 
implies that the new capabilities of the VLA and ALMA 
offer a unique opportunity to detect a much larger pop¬ 
ulation of radio sources in young stellar clusters. New 
observations with improved sensitivity and better angu¬ 
lar resolution will provide crucial information about the 
origin and nature of the radio emission, and they will 


reveal how radio and X-ray phenomena are connected. 

Furthermore, the presence of multiple variable radio 
sources would have important implications for interfero¬ 
metric imaging, since the classical techniques assume a 
constant sky. A more accurate determination of the radio 
flare rate would help to understand how this variability 
can affect the upcoming VLA and ALMA observations 
in young stellar clusters. 
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